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ABSTRACT 

We report the detection of luminous CO (3-2) line emission in the halo of the z=2.6 
radio galaxy (HzRG) TXS0828+193, which has no detected counterpart at optical to 
mid-infrared wavelengths implying a stellar mass < fewxlO 9 M and relatively low 
star- formation rates. With the IRAM PdBI we find two CO emission line components 
at the same position at ~80 kpc distance from the HzRG along the axis of the radio 
jet, with different blueshifts of few 100 km s _1 relative to the HzRG and a total 
luminosity of ~ 2x 10 10 K km s _1 pc 2 detected at a total significance of ~ 8<j. HzRGs 
have significant galaxy overdensities and extended halos of metal-enriched gas often 
with embedded clouds or filaments of denser material, and likely trace very massive 
dark-matter halos. The CO emission may be associated with a gas-rich, low-mass 
satellite galaxy with very little on-going star formation, in contrast to all previous CO 
detections of galaxies at similar redshifts. Alternatively, the CO may be related to a 
gas cloud or filament and perhaps jet-induced gas cooling in the outer halo, somewhat 
in analogy with extended CO emission found in low-redshift galaxy clusters. 
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1 INTRODUCTION 

The most vigorous starbursts in the Universe occurred in 
massive galaxies during the most active phase of galaxy evo- 
lution and AGN activity, at redshifts z~ 2—3. These galaxies 
formed most of their stellar mass of a few xlO 10-11 Mq in 
short bursts of few xlOO M yrs with star-formation rates of 
several 100 yr" 1 (e.g.. ISmail et aDl2002l : lArchibald et al.1 
l200ll : lReuland et al.ll2004h . Luminous CO emission observed 
at millimeter wavelengths is typically viewed as the most 
direct sign of the immense reservoirs o f cold molecular g as 
necessary to fuel these starbursts (e.g. jGreve et al.ll2005h . 

Powerful high-redshift radio galaxies (HzRGs) may 
host the most extre me starbursts at high redshift (e.g., 
ISevmour et al. 2008), seen in a short, but critical phase 
of their evolution dominated by strong AGN feedback 
(jNesvadba et alJfeOOd I2007L l2QQ8h . Bright K-band magni- 
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tudes suggest HzRGs are a mong the most massive galaxies 
at all cosmic epochs (e.g., iDe Breuck et al.l l2QQ2h , a con- 
clusion recen tly confirmed th rough rest-frame near-infrared 
photometry. ISevmour et~ai1 ([2007T ) find stellar masses of 
several 10 11 M©, factors of a few larger than typical 
mass es of submillimete r galaxies at similar redshifts (10 10 ' 5 
Mm. ISmail et al.l [2004 V Lumin ous CO emission has been 
found in several HzRGs (e.g., Papado poulos et all [2000; 



IDe Breuck et al.ll2005l : lKlamer et alJ2005h. and rec ently also 
in a satellite galaxy of a HzRG (|lvison et al.ll2008h . 

HzRGs reside in particularly rich environments, and 
are o ften surrounded by several 10s of companion galaxies 
(e.g. iLe Fevre etHlll996l : iKurk et al.ll2004l : IVenemans et al.1 
2007) as well as extended halos of ionized and neutral gas. 
Faint, diffuse Lya emission extends to radii well beyond the 
inne r halo, where the g a s is strongly disturbed by the radio 
iet. lVillar-Martm et al.l ([20021 120031 ) trace ionized gas out to 
radii of > 100 kpc with relatively quiescent kinematics and 
CIV emission line ratios implying near-solar metallicites out 
to large radii. Deep Lya absorbtion troughs reveal neutra l 
gas, likely in clouds or filaments (e.g.. Ivan Oiik et al.|[l997h . 
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Figure 1. CO(3-2) line image, integrated over the full wavelength 
range of the two components. Line emission is not spatially re- 
solved with a h" beam. Thick yellow and thin white contours show 
the non-thermal continuum observed at 3mm with the PdBI and 
at 8.4 GHz with the VLA, respectively. 



Using the IRAM Plateau de Bure Interferometer we de- 
tected luminous CO (3-2) line emission in the halo of the 
z=2.6 HzRG TXS0828+193, with a luminosity of 2 x 10 10 
K km s _1 pc 2 . Deep photometry from the rest-frame UV to 
mid-infrared including MIPS 24/xm imaging does not reveal 
a counterpart within the 5" beam, implying a very small as- 
sociated stellar mass and low star- format ion rates. These are 
very unusual properties for a high-redshift CO emitter, and 
we discuss possible scenarios for its nature. Throughout the 
paper, we adopt a flat Ho =70 km s _1 Mpc -3 concordance 
cosmology with Qa — 0.7 and — 0.3. 



2 OBSERVATIONS AND ANCILLARY DATA 

We observed TXS082 8+193 with the IRAM P lateau de Bure 
Interferometer PdBI (jGuilloteau et al.l ll992) in the D con- 
figuration. At z=2.6, CO(3-2) falls at 96.6 GHz and into the 

3 mm atmospheric window. We reached an 0.3 mJy/beam 
rms and a beam size of 5.3" x 4.6" (42 x 37 kpc at z= 2.6). 
On-source integration time was 12.9 hrs under normal to 
good conditions with 6 antennae and system temperatures 
< 150 K. Data were calibrated using the CLIC package and 
with MWC349 as flux calibrator. We combined both polar- 
izations and rebinned the data to a resolution of 30 km s _1 . 
The PdBI receivers covered a window of ~ 2600 km s _1 . 

We also include deep Palomar WIRC K-band imag- 
ing, and SPITZER IRAC 3.6/xm and MIPS 24/rni pho- 
tometry, as well as archival HST WFPC2 F606W imag- 
ing. C. Carilli kindly provided his VLA 8.4 GHz A-arra y 
map of TXS0828+193 first published in lCarilli et all (jl997f ). 
We detected the millimeter continuum of the radio core of 
TXS0828+193. Thus, we could align the PdBI data with 
our line- free K band continuum image of TXS0828+193 
([Nesvadba et a l. 2008), assuming that the AGN is at the 
center of the galaxy. The K band data then served as a ref- 
erence to align all other images. 



3 CO EMISSION IN THE HALO OF A HZRG 

We identify luminous CO (3-2) emission in the halo of 
TXS0828+193 at a distance of - 10" (80 kpc) south-west 
from the radio galaxy (Fig. [T]). The 5" beam corresponds to 
an upper size limit of ~ 40 kpc. The emission line region is 
aligned with the axis of the radio jet but 2.5" (20 kpc) WSW 
in projection from the radio lobe, and hence at a larger ra- 
dius from the central galaxy. Integrating over the full line 
we reach 8a significance. To further ensure the robustness 
of our measurement, we split the data into two subsamples 
with consistent results in each subset. 

The integrated spectrum is shown in Fig. [2] 
We find 2 compact components, labeled TXS0828+193 
SW1 and SW2, respectively. They have the same 
spatial position, but different blueshifts of Avswi=- 
200=b40 km s" 1 and Av S W2=-920±70 km s" 1 rela- 
tive to the HzRG, respectively. The systemic velocity 
was estimated from rest-fra me optical integ r al-fiel d spec- 
troscopy of TXS0828+193 ([Nesvadba et al.1 l2008h . SW1 
and SW2 have line widths FWHM S wi=310±270 km 
s" 1 and FWHM SW 2=340±270 km s~\ respectively. In- 
tegrated fluxes are Ico,swi=0.23=b0.08 Jy km s _1 and 
Ico,sw2=0.24±0.06 Jy km s _1 , respectively, and correspond 
to luminosities of C' co —9x 10 9 K km s _1 pc 2 per component. 
(We assumed C / 3 _ 2 —^i-o and 7*32=1) . 

Molecular gas mass estimates of high-redshift galaxies 
depend on the assumption that the conversion factors from 
CO to H2 ( "X factors") established at low redshift will apply. 
The X-factor of Xu = 0.8 M©/(K km s _1 pc 2 ) appropriate 
for ULIR Gs yields estimates — 5x lo wer than the Milky Way 
X-factor dDownes fe SolomoiJ[l998h . With Xu, the I S wi ~ 
Ssw2 ~ 0.25 Jy km s _1 per component correspond to 7x 10 9 
M in cold gas per component, or 1.4 x 10 10 M in total. 

We will now discuss two possible scenarios for the na- 
ture of the CO emission in the halo of TXS0828+193, 
namely, that it may be associated with a satellite galaxy, 
or that it may be related to gas clouds or filaments within 
the gas-rich halo of TXS0828+193. 



3.1 An extremely gas-rich satellite galaxy? 

We searched for the stellar continuum of putative galaxies 
associated with SW1/2 in our set of images with rest-frame 
wavelengths between ~ 1700A and 7/im (J5J). SW1/2 was 
not detected in any of the data sets (Fig. [3}. This is in strong 
contrast to the companion of the z=3.8 HzRG 4 C60.07, that 
was d etected in all bands in a similar study of llvison et al.l 
(2008|). We use the deep K-band photometry with a 3a limit 
of K3 a =23.7 ma p in a \" aperture and th e population syn- 
thesis models of Bruzual & Charlotl £2003), to place an up- 
per limit on the stellar mass. Continuous star- format ion his- 
tories with ages between 5 x 10 7 yrs and 2 x 10 9 yrs (implying 
a formation at z<10), and extinctions Ky = 1 — 5 mag cor- 
respond to a maximum of - 3x 1O 9 M in stellar mass. This 
covers more than the range of extinctions found for dusty 
submillimeter galaxies at 2 and low-redshift ULIRGs 
typically associated with strong CO emission (Ay < 2 mag 
ISmail et al.ll2004l : IScoville et alfeoooh . These extinctions are 
derived from the integrated photometry of the galaxies, as 
appropriate for our purposes. Extinctions along individual 
sight lines and into a starburst may be significantly higher. 
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Figure 2. Integrated spectrum of the CO emitter south-west 
from TXS082+193. The yellow histogram marks the data, with 
a Gaussian fit of two components overlaid (black line). The 
blueshifted and redshifted components are detected at 3cr and 
4cr significance, respectively. 
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Figure 4. Lyo; spectrum of TXS0828+193, graciously provided 
by M. Villar-Martin with the CO(3-2) components SW1 and SW2 
shown as red and green ellipses, respectively. Positions and sizes 
correspond to the distance from the HzRG and beam size, and 
to the relative redshift and FWHM line width of the two compo- 
nents, respectively. 



The r.m.s. of ~ 100 /iJy i n our MIPS 2 4/im image is well 
below the fluxes measured bv lPope et al. (2008) for submil- 
limeter galaxies at similar redshifts, which are in the range 
200-500/iJy. This allows us to set constraints on the star 
formation, because the filter covers the 6.2/xm PAH band, 
and more than half of the 7.7/xm band at the redshift of 
SW1/2. With the MIPS non-detection, it appears unlikely 
that SW1/2 is forming stars at the prodigeous rates of sev- 
eral 100 M typically observed in SMGs. 

We will now estimate a dynamical mass for SW1/2. 
With a beam size of 40 kpc, we do not know whether 
the CO line emission may be associated with one or with 
two galaxies. If SW1/2 represent a double- horned profile 
of a single, roughly virialized, rotating galaxy (as often 
assumed for SMGs), we can estimate a mass by setting 
Md yn = (v/ sini) 2 Rik pc /G with circular velocity v =350 
km s _1 , inclination z, radius Rikpc, a nd gravitation a l con - 
stant, G. The radius R is given in kpc. lTacconi et all ([2008) 
use radii ^2—5 kpc for their mass estimates, which would 
imply Mdyn ~ 0.6 — 1.5 x 10 11 M© for an edge-on disc, 
and, perhaps more realistically, 2-3 x higher masses for a 
more average inclination. Thus, SW1/2 would have a mass 
in the typical range of submillimeter galaxies or powerful ra- 
dio galaxies, which are 2 mag or brighter (|Smail et alj 
120041 : iDe Breuck et alJl2002h . This would also significantly 
exceed the baryonic (gas and stellar) mass of few xlO 10 M . 
If alternatively, we assume that SW1 and SW2 are associ- 
ated with two different galaxies in the halo of TXS0828+193 
(a plausible assumptio n given the 40 k pc covered by the 
beam), then, following iNeri et all (|2003h , we estimate a dy- 
namical mass of Mdyn = 4 x 10 9 Rikpc M© per galaxy for 
a FWHM=300 km s _1 line width of each component. As- 
suming a radius of a few kpc, the dynamical mass estimate 
will be lower than the molecular gas mass by factors of a 
few, except if we assume that both galaxies are seen within 
a few degrees from being face-on, which does not appear 
very likely. Each of these estimates relies on the assumption 
that the gas is approximately virialized. This is a common 
assumption in CO emission line studies at high redshift, but 
whether it is justified has yet to be proven. In fc|3.2l we dis- 
cuss a sc enario where the vi rial assumption would not apply. 
Likewise llvison et al.1 ([2008) raised doubts as to whether this 



assumption is always justified in the context of high-redshift 
galaxies. 

In these "galaxy" scenarios, we also need a mechanism 
to suppress star formation in the cold gas traced by the CO. 
If the CO line emission arises from a disc with a few kpc in 
radius, then the observed gas mass of 1.4 x 10 10 M corre- 
sponds to a surface mass density of few x 1000 Mp pc -1 . 
Following the Schmidt-Kennicutt relation (|Kennicuttl ll998) 
between gas surface density and star formation intensity, 
SFI, we expect SFI = few x 10 M yr _1 kpc -2 . Averaging 
over the size of the disc, this corresponds to a total star- 
formation rate of several 100 Mq yr -1 . This is in the typical 
range of submillimeter galaxies, but in contradiction with 
our non-detection at 24/xm. Likewise, with star- format ion 
rates of fewxlOO M© a stellar mass of fewxl0 9 M© would 
be built in fewxlO 7 yrs, so that SW1/2 would have to be in 
a very special, very young stage of the starburst if it was a 
galaxy. 

IPapadopoulos et al.l (|2008h recently found luminous, 
but excited CO line emission in a nearby radio galaxy, 
3C293, which does not seem associated with a starburst. The 
same is suggested by Spitzer observations of a small number 
of nearby galaxies with strongly enh a nced, mid-infra r ed Hg 
line emission (e.g., iQgle et al.1 120071 ). iGuillard etHI (|2009h 
argue that this gas may be heated through the dissipation of 
kinetic energy. In these cases, it is likely that the energy was 
injected by an external mechanism (a merger or AGN). We 
will in the following propose a somewhat related scenario, 
where the nearby radio lobe may have induced the collapse 
of gas within the halo of TXS0828+193. 



3.2 Cold gas in the halo? 

TXS0828+ 193 was the first HzRG where an outer halo 
was found ([ Villar- Mart m et al.l l2002). which extends beyond 
the turbulent, luminous inner emission line region that is 
likely powered by energy relea s ed from the powerful AGN 
(e.g., Villar -Martm et al. 1999: iNesvadba et al.l [2008). Dif- 
fuse gas in the outer halo is fainter, and has more mod- 
erate velocity gradients and line widths, which may indi- 
cate r otation in the potential of the u nderlying dark- matter 



halo (I Villar- Martin et al] 



20031. l2006h, or perh a ps ga s in- 



fall ( Humphrey et al.ll200/t ). IVillar-Martm et al.1 ([2002) find 
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Figure 3. (left to right:) HST WFPC606W, Palomar K band, Spitzer IRAC 3.6/im and MIPS 24/im photometry of TXS0828+193. 
Thick blue contours show the position of the SW1/SW2, thin red contours mark radio jets. SW1/SW2 is undetected in all bands. 



CIV emission in the outer halo of TXS0828+193 and sus- 
pect that metallicities may be up to nearly solar, perhaps 
representing gas that has previously been driven out from 
the central galaxy. 

Fig. S] shows the relative velocity and radial distance 
of SW1 and SW2 from TXS0828+193 relative to the two- 
dimen sional Lya longslit spectrum of IVillar-Martm et al.l 
(2002). The slit of the Lya spectrum falls near the CO po- 
sition, but since the slit width, 1", is significantly smaller 
than the beam of 5", we cannot infer whether it covers the 
position of the CO emitter. Lya emission is traced almost 
out to the 24 kpc radial distance of SW1/2 from the HzRG. 

SW1 and SW2 are blueshifted relative to the ra- 
dio galaxy, and are also on the side of the radio galaxy 
where the diffuse Lya halo is blueshifted (the luminous in- 
ner halo is redshifted). The kinematic properties of SW1 
(FWHM=310±270 km s~\ z=2.5761±0.0008) are very sim- 
ilar to those of the diffuse ionized gas (FWHM=590±60 km 
s _ \ z=2.5741±0.0008), with a relative blueshift of SW1 
relative to the Lya emission of 170 ±60 km s" 1 . SW2 
(FWHM= 340 ± 270, z=2.5678±0.0007) is blueshifted by 
890 km s _1 relative to the diffuse Lya emission, but has 
a similar line width within large uncertainties. Given the 
low stellar mass possibly associated with the CO line emis- 
sion (§ 13. ip , this may suggest that the CO emission orig- 
inates from clouds or filaments in the diffuse halo gas of 
TXS0828+193, which has metallicities of up to ab out solar 
(|Villar-Martm et al.ll2002h . lDe Breuck et al.l (|2QQ3h detected 
CO emission at the redshift of a Lya absorber near the z=3.1 
HzRG B2 2330+3927, and proposed a similar scenario, but 
did not have the spatial resolution to directly measure posi- 
tional offsets between the radio galaxy and CO emission. 

Several mechanisms may plausibly influence the halo 
gas including minor or major mergers, or powerful outflows 
from starbursts and AGN. Each of these mechanisms may 
sweep up and accelerate halo gas over timescales of a few 
x 10 8 yrs similar to those suggested by the relative velocity of 
the CO and distance to the radio galaxy. The close proximity 
of the emitters to the radio hot spot and alignment with the 
jet axis is however suspicious. The emitters are very close 
to the jet axis, and within a projected area of ~300 kpc 2 
from the hot spot, whereas our data have a half-power beam 
width covering a total of 135,000 kpc 2 . If this is not due to 
mere projection effects, then weak shocks produced by the 
expanding radio source may play a role in compressing and 
exciting the gas (the radio hot spot is only about 2.5" or 20 



kpc away, and we may not detect the low surface brightness 
radio plasma). In turn, the interaction with dense gas may 
be enhancing or even triggering the radio hot spot. 

Extended cold molecular gas is found in some mas- 
sive cooling-flow cluste rs at low redshift (e.g., lEdgej l200ll ; 
ISalome fe Com bes 2004), where molecular gas forms along 
the edges of X-ray cavities inflated by the radio jet and out- 
side the volume filled by the radio plasma. Similarly, CO 
line emission in the halos of HzRGs may trace the edges of 
cavities inflated by the radio lobes. Within the large uncer- 
tainties, the -300 km s" 1 FWHM of SW1/2 are not very dif- 
ferent f rom the line widths in the diffuse CO in the Perseus 
cluster (jSalome et alJl2008h . 

We can use the observed surface brightness of the faint 
Lya emission in the halo of TXS0828+193 and the observed 
CO luminosity to investigate whether TXS0828+193 falls 
near the correlation between the luminosity of the molecu- 
lar and ionized gas found in local cooling-flow clusters (Edge 
l200lh . Adopting a flux conversion Cl V (x — 13 x Cho. be- 
tween Lya and Ha luminosity (where we neglect extinc- 
tion), we find a strict upper limit on the Lya luminosity of 

£ L = 1.4 x 10 42 erg s" 1 cm -2 within the 20D" ar ea of 

I 1 1 1 

the beam. Translating the CO gas mass of Edge (2001) into 
a CO luminosity, we find that the halo of TXS0828+193 
falls only factors of a few below the expected value found 
in local cooling- flow clusters. Allowing for different physical 
conditions, gas distributions, and AGN properties (HzRGs 
host powerful AGN), we may well be seeing a fundamentally 
similar phenomenon. 

However, the brightness temperatures and spatial dis- 
tribution of the extended CO emission in low-redshift clus- 
ters are significantly different. Low brightness temperatures 
suggest low gas filling factors, and much of the gas is concen- 
trated towards the central radio galaxy. This may arise from 
different properties of the diffuse cluster gas at high and low 
redshift, which is cold for HzRGs (with embedded filaments 
or clouds of few xl O 9 M and more in dense neutral gas 
Ivan Oiik et al.l Il997h and hot and rarefied in low-redshift 
clusters. In fact, in HzRGs at z^ 2-3 we may be witness- 
ing the processes that led to the rapid heating and entropy 
enhancement of cluster gas through AGN feedback, which 
seem necessary to explain the temp erature profiles of mas- 
sive X-ray clusters at low re dshift (|Nath &; Rovchowdhurvl 
120021 : iMcCarthv et alJ l2008h . Interestingly, Lya absorbers 
are only found in the halos of HzRGs with relatively small, 
and likely rather young, radio sources. This suggests that 
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large amounts of dense gas (and perhaps dust) may be 
present in the halo at radii that are not yet affected by 
mechanical heating from the radio source. The approaching 
radio jet of TXS0828+193 may have contributed to trigger- 
ing the collapse and forming SW1/2. Similar processes may 
ultimately lead to positive AGN feedback and jet-triggered 
star formation in some cases, if the cold molecular gas will 
relax and form stars over sufficiently short time scales. 



4 SUMMARY AND CONCLUSIONS 

We discuss the nature of luminous CO (3-2) line emission in 
the halo of the radio galaxy TXS0828+193 at z=2.6. The CO 
emission resembles that of submillimeter galaxies, but we do 
not detect continuum emission from SW1/2, including 24/xm 
MIPS imaging, which covers the PAH bands at z=2.6. For a 
gas disc in a galaxy we would expect strong star formation 
if the Schmidt-Kennicutt law roughly applies. 

Alternatively, SW1/2 may represent a cloud or filament 
in the halo, maybe related to neutral, dense Lya absorbers 
observed near some HzRGs. The approaching radio jet of 
TXS0828+193 may have contributed to triggering the col- 
lapse and exciting the gas. This is somewhat in analogy with 
diffuse CO emission in low-redshift clusters, but the ambi- 
ent gas properties will likely be very different at z=2.6. In 
either case, SW1/2 does not appear to be an 'ordinary' high- 
redshift CO emitter, and further observations are necessary 
to differentiate between the two scenarios. This suggests that 
CO observations of the high-redshift Universe with the re- 
furbished PdBI and soon with ALMA, will reveal a rich, and 
multi-faceted picture of the early Universe. 
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